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Abstract 

We studied possible evolution of rotational period and magnetic field 
of the X-ray source RX J0720. 4-3125 assuming this source to be an iso- 
lated neutron star accreting from interstellar medium. Magnetic field of 
the source is estimated to be 10 6 — 10 9 Gs (most probably ~ 2 • 10 8 Gs), 
and it is difficult to explain the observable rotational period 8.38 s with- 
out invoking hypothesis of the magnetic field decay. We used the model 
of ohmic decay of crustal magnetic field. The estimates of the accretion 
rate (1CT 14 - lO _16 M /yr), velocity of the source relative to the interstellar 
medium (10 — 50 km/s), the neutron star age (2- 10 9 — 10 10 yrs) are obtained. 



(remark added in proof: As we found out in July 1997 when this paper was 
already submitted, common results were recieved independently by John 
C.L. Wang // Astrophys. J., 1997, V.486. P.L119) 
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1 Introduction 



Recently isolated neutron stars (INSs), which are not observed as radiopul- 
sars, became of great interest both for observers and theorists. The idea 
of observations of such sources was proposed more than 25 years ago [6], 
and in 1991 Treves and Colpi [8] suggested, that INSs, accreting from inter- 
stellar medium (ISM), can be observed in great number in UV and X-rays 
by ROSAT. Here we present the results of the work on the source RX 
J0720.4-3125 observed by Haberl et al. [12]. In [12] it was proposed, that 
RX J0720. 4-3125 is an accreting INS with the rotational period 8.38 s. 

There are four possible stages of an INS in low density plasma: Ejector 
(E), Propeller (P), Accretor (A) and Georotator (G). The stage is deter- 
mined by relations between specific radii: R[- light cylinder radius, R s t~ 
stopping radius, Rq = { 2 % M )~ the radius of the gravitational capture and 

CO 

Rco = (^rrO 1 / 3- the corotation radius. 

As a result we have two critical periods: Pe and Pa, separating different 
stages. If p < Pe, we have an Ejector (i.e. a pulsar), if Pe < p < Pa, 
the NS is on the Propeller stage, and if p > Pa and R s t < Rg - we have 
an accreting NS. It is possible, that p > Pa, but R s t > Rg- In this case a 
geo-like magnitosphere is formed, and we call it Georotator. 

For describtion of period evolution it is useful to use a gravimagnetic 
parameter, y, (it is described in [4]). 

On the figure 1 three examples of evolutionary tracks of an INS on p-y- 
diagram are presented. All of them are ended at the Accretor stage. 

I. E — ► P — ► A- evolution of an INS in the ISM with the constant density 
without magnetic field decay. 

II. E — ► P — ► A — ► P — > A- evolution of an INS, passing through a 
giant molecular cloud without magnetic field decay. 

III. Evolution with magnetic field decay. 

We are especially interested in the third case. 

2 Analitical estimates of the X-ray pulsar's pa- 
rameters 

In this part we show mainly estimates which were not included into the 
article [3]. 

Lets consider the situation, when the Alfven radius, Ra, and the coro- 
tation radius, R co , are equal. This equation gives us the accretion period, 
Pa: 
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p,-6.io^ 7 p:^ 6 7 (^)" 1/7 s (i) 

here - magnetic momentum in units 10 30 Gs • cm 2 , /?_24 _ density of the 
ISM in units 10~ 24 g/cm 3 , Voo 6 - velocity of the INS relative to the ISM in 
units 10 6 cm/s. 

During accretion p > Pa- If these two periods are equal: 

6/7 _ 8.38 3/7 9/y (M_\ 1117 

So, /i3o = O.OOTp^^oo^ 2 ^ for = p = 8.38 s (the period of 

RX J0720.4-3125). For the NS's radius R = 10 km we have B < 7 • 10 9 Gs. 

If we use the hypothesis of the acceleration of an INS from the turbu- 
lizated ISM, we have an equation, which gives an estimate of the magnetic 
field of the INS based on the equilibrium period, P eq , [3, 5]: 

P eq = 2355^/3^34/3^2/3 ^3/3,-2/3 ( g (2) 







Here I45 - the momentum of inertia in units 10 45 g • cm 2 , vt 6 — the tur- 
bulent velocity in units 10 6 cm/s. 
So we have: 



M30 



8.38 N 3 / 2 1/2 1/2 _ /2 / M 
2355 J *« 745 ^ 24 ^ ^ I Mo 



And finaly 5 « 2.1 • 10 8 Gs. 

In [3] we showed, that the magnetic field can't be less than 10 5 — 10 6 Gs, 
because in the opposite case we can't observe pulsations of the X-ray flux. 

Lets suppose, that the INS was born with a period about 0.01 - 0.02 s 
(the exact value is not very important, because the only necessary condition 
is Pinitial « 8.38 s, and initially the INS was on the Ejector stage). So, the 
star had to decelerate till 8.38 s. Lets estimate the time of this deceleration. 

The time of deceleration is determined by the final, not by the initial 
period! 



p E ~ nktf/^pztfv^s 



dloj [i 2 

-h-yr, Ri = c/cu, lo = 2ir/p 



dt l Rf 
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Figure 1: P-y diagram 
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Figure 2: Changes of the surface magnetic field of an isolated neutron star 
with time for the model of standard cooling. Curves 1, 2, and 3 correspond 
to the initial depths of the current layer, 10 11 , 10 12 , and 10 13 g cm~ 3 , 
respectively. The solid curves correspond to Q = 0.001; the dashed curves, 
to Q = 0.01; the dot-dashed curves, fe Q = 0.1. 




Figure 3: The evolutionary tracks of the neutron star for the accretion 
rates M = Hr^-Moyr -1 (a) and M = 1(T 16 M yr" 1 (b). The model 
parameters are described in the text. The dashed lines correspond to p = 
Pe; the dot-dashed lines, to p = Pa- The dashed line in Fig. 2a shows 
for the second track the neutron star evolution with no acceleration in the 
turbulized intestellar medium. The numbers near the marks in tracks denote 
the logarithm of the neutron star age in years. The observed radio pulsars 
are indicated by dots. 
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So for Ap = p: 



l 7 Ap = k /{^f 

p2 c 3p3 

Ip 2 c 3 

At - 



= 3 • W 7 P E K 1 ^30 h 5 yrs = 
= 3 • W 9 k 1 T l/2 n^pZ l 2 / 4 2 v^ l f .h 5 yrs 

For ^30 < 0.01 At > 3 • 10 n yr,s >> t Hubb i e . I.e. the star couldn't 
initially have low magnetic field. 

As far as the NS is on the Accretor stage, its period should be much 
greater than the ejection period, Pe- 

P > Ppropeller > Pe 

Pe ~ 10s for standart parameters. It means, that if the star is on the 
accretor stage with p=8.38 s, the field is much less than the standart value: 
B « 10 12 Gs. 

Lets try to estimate the characteristic time of acceleration and deceler- 
ation of such a NS: 

Ilo 



20yrs ■ hsv^p 

-24 



Mv t R G 



5 „-l / P \ 1 / V t 



°°^- 24 Vl0 5 s; \W 6 cm/ S/ 

where t su ,t s d - characteristic times of the period changes. 

For p=10 s t su = 2- 10 5 yrs. If v t is less than 10 6 cm/s, than t su is greater. 

8.38 
210 b -310 Y 



So, p » p/t su orMlw ~ ^- 12 s/s 



3 Calculations of the magnetic field decay 

So, we can say, that the magnetic field of the source dissipated, and the 
characteristic time of the dissipation , td, was short enough: td < tE (tE~ 
Ejector's time). Because of that we have rapidly rotating INS on the Accre- 
tor stage. 

The detailes of calculations one can find in [3]. 
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We calculated the magneto-rotational evolution of the star with the 
mass M = 1.4M Q for accretion rates lO -15 M0/yr and 10 _16 MQ/yr passing 
through the ISM with the density p = 10~ 24 g/cm 3 . These rates correspond 
to velocities ~ 20 and ~ 40 km/s. 

The magnetic field decay was studied by different authors (see [10]). For 
conductivity we used formulae from [13] and [2]. Initially the magnetic field 
was localized in the surface layer. Low accretion rates don't have much 
influence on the magnetic field decay [1, 9]. 

On the figure 2 we show the decrease of the surface magnetic field for 
different parameters. We used the model of the NS structure from [11]. 

On figure 3 the evolutionary tracks for the accretion rates lO _15 M0/yr 
(fig. 3a) and 10 _16 M©/yr (fig. 3b) are shown. 

4 Conclusions 

Observations of the accreting INSs can be an important test for theories of 
the magnetic field decay. In that case we have a very " pure" example of the 
decay, because a lot of effects of the interaction of the magnetic field, the 
surface of the NS and surrounding plasma are not important. 

Observations of the source RX J0720. 4-3125 [12] showed the existence 
of the INS with the rotation period which can be easily explained with the 
assumption of the magnetic field decay. It can be important for estimates 
of the total number of observable INS and for their appearence as periodic 
X-ray sources [7]. 

We want to thank dr.F. Haberl for the information about the source, 
drs. M.E. Prokhorov, V.A. Urpin, D.G. Yakovlev and V.M. Lipunov for 
helpful discussion, and participantes of the summer school "The many faces 
of neutron stars" , where the idea of this article was born. 
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5 References 

1. Zdunik J.L., Haensel P., Paczyhski B., Miralda-Escude J.// Astrophys. 
J. 1992. V.384. P.129. 

2. Itoh N., Hayashi H., Kohyama Y. // Astrophys. J. 1993. V.418. P. 405. 

3. Konenkov D.Yu, Popov S.B.// Astronomy Letters. 1997. V.23, P. 200 
(|astro-ph/9707318|) 

4. Lipunov V.M./ / Astrophysics of neutron stars, Springer, 1992. 



8 



5. L ipunov V.M., Pop ov S.B.// Astron. Zhur.. 1995. V.72, P. 711. (see 
also |astro-ph/9609185| ) 

6. Ostriker J.P., Rees M. J., Silk J. // Astrophys. J. Letters 1970. V.6, 
P. 179 

7. Popov S.B. //Astron. Circ 1994. No.1556. P.l. 

8. Treves A., Colpi M.// Astron. and Astrophys., 1991, V. 241, P. 107. 

9. Uprin V., Geppert U., Konenkov D.// Astron. and Astrophys. 1996. 
V.307. P.807. 

10. Urpin V.A., Muslimov A.G.// Astron. Zhur. 1992. V.69. P.1028. 

11. Friedman B., Pandharipande V.R.// Nucl. Phys. 1981. A361. P.502. 

12. Haberl F., Pietsch W., Motch C, Buckley D.A.H.// IAU Circ. 1996. 
No. 6445. 

13. Yakovlev D.G., Urpin V.A.// Astron. Zhur. 1980. V.24. P.303. 



9 



